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ABSTRACT 


The purpose of this thesis is to develop a means by which 
the transient response of a servomechanism may be rapidly 
evaluated. The approach consists of developing a method of 
rapidly determining the value of several transient response 
parameters so that the system transient response may be sketched. 
The parameters selected for this purpose are time and magnitude 
of peak overshoot, rise time and settling time. 


To permit evaluation of these parameters systems are 
classified according to the number and type of singularities 
in their closed loop performance functions. Systems considered 
are underdamped systems with two or three singularities. For 
each class of system, loci of constant values of overshoot and 
constant values of the three time parameters expressed in terms 
of an artificial time are mapped onto a complex plane. Para- 
meter loci for systems with only two singularities are mapped 
directly onto the s-plane. Loci for systems with three singu- 
larities are mapped onto an auxiliary complex plane. 


The values of the four transient response parameters for 
any system of a particular class can be read from the families 
of loci for that class. The magnitude of overshoot can be 
used directly and the three time parameters can be easily con- 
verted into real time. The initial slope of the transient 
response and the slope at the time of peak overshoot are either 
known or easily determined. Also the transient oscillating 
frequency is obtainable directly from the system closed loop 
performance function. Using this information the transient 
response of a particular system can be sketched with considerable 
accuracy. This sketch can then be used as a basis for rapid 
System evaluation. 
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Chapter I 
INTRODUCTION 


In the past two decades the importance of servomechanisms 
and automatic control systems in general has increased tre- 
mendously. Accordingly, much effort has been expended in de- 
veloping methods of analysis and design. The result is an 
abundance of ideas and techniques available to the designer 
today. Most of these can be found in any one of a number of 
current automatic control system textbooks. 

Basically there are two approaches to analysis and design, 
the frequency response approach and the time domain approach. 
Either of these techniques may be used with some facility in 
designing systems to meet both time domain and frequency re- 
sponse specifications. However, if time domain specifications 
are to be met, a time domain analysis must eventually be made. 
This is true regardless of which technique is used in the 
actual design. 

The usual time domain specification is one involving 
transient response, that is, tne time response of the system 
to a suddenly applied input signal. The most common input 
signal used is the unit step. Hereafter we will use the term 
transient response exclusively to mean the time response of a 
System to a unit step input signal with all initial conditions 
equal to zero. 
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The problem of obtaining the transient response of servo- 
mechanisms is one of considerable magnitude. The transient 
response of linear second order systems has been studied ex- 
tensively and is well known. The transient responses of high- 
er order systems and second order systems that have been 
complicated by addition of a zero are more difficult to obtain 
and are 1ess well known. The exact transient response of any 
linear system can always be obtained by factoring the charac- 
teristic equation, solving for the output as a function of 
time and plotting this function. However, this procedure in- 
volves considerable labor and is quite time consuming. 

One method of attempting to avoid such a calculation is to 
consider the complex roots of the characteristic equation to 
be dominant and use a second order approximation of the more 
complicated system. As the value of the real roots and zeros 
become smaller with respect to the undamped natural frequency 
of the system this approximation ceases to be valid. When the 
second order approximation ceases to be valid another approxi- 
mation developed by Chu (1) may be employed. This method re- 
duces the labor involved in obtaining the transient response 
by applying corrections to the second order equations for 
zeros and other roots of a system. Unfortunately this approxi- 
mation also becomes invalid as the values of the real roots 
and zeros become even smaller. 

If the frequency response approach is used, obtaining the 
transient response is usually more difficult because the system 
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Characteristic 


equation is not known in factored form. In 


this case Sraphical integration Schemes such as one developed 


by Floyd (2) are sometimes used to obtain an approximate 


transient response. These methods also require considerable 


labor and time, 


In this thesis we restrict the discussion to linear 


Systems which exhibit the property of unity feedback, that is, 


Systems for which the steady state Output becomes unity when 


a unit step input signal is applied, Notice that absolute 


Stability is implied, 


With this restriction we can make some genera] obser- 


vations about the transient response of any such System, 


a. 


At zero time the value of the output 
1s zero, 


At a Sufficiently large time the 
Value of the output becomes unity, 


For critically damped and ovepr- 
damped systems the output increases 
monotonically with time until ac 
reaches its Steady state Value of 
unity. This case is of little 
practical interest, 


A more common and more practical 
case is that of an underdamped 
System, In this case the output 
rises to a value greater than 
unity, and may or may not oscillate 
about unity before assuming its 
Steady state Value, 
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Based on the above assumptions and observations we may 


define the following transient response parameters: 


Peak Overshoot - M The amount by which the maxi- 
eee ee 


mum value of the output exceeds the steady state value. 


Time of Peak Overshoot - t The time at which the 
i. a ee 


maximum value of output occurs. 


. Rise Time - ty The time at which the value of the out- 


put is unity for the first time. 


Settling Time - t. The mathematical expression for the 


transient response of systems which have only one pair 
of complex roots may consist of three parts. These are 
a constant equal to the steady state value of output, 
decaying exponential terms, and a damped sinusoidal 
term. The settling time is the time at which the sum 
of the error due to the exponential terms plus the 
height of the envelope of the damped sinusoid is equal 


to two per cent of the steady state output value. 


A sketch of a typical transient response with the above- 


defined parameters indicated is shown in Figure I-1l. 


If the four transient response parameters defined above 


and the transient oscillating frequency of a system are known, 


a Sketch which very closely approximates the actual transient 


response can be constructed. A means of rapidly determining 
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FIGURE I-1 
TYPICAL TRANSIENT RESPONSE 
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these parameters would then be a valuable aid to the de- 
Signer in evaluating the transient performance of a system 
without computing the exact transient response. 

In this thesis a technique which permits rapid graphical 
determination of these parameters by inspection is developed 
for certain classes of systems. The technique involves the 
mapping of lines of constant value of these parameters onto 
a complex plane for a particular class of system. The values 
of these parameters for all systems of that class can then 


be read directly from the complex plane plots. 
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Chapter II 
BASIS OF THE MAPPING TECHNIQUE 


The determination of the transient response parameters of 
a given system is most easily approached through the per- 
formance function, the singularities of which can be displayed 
om a complex plane. In this approach the Laplace transform 
with the associated complex function theory is used as a 
mathematical tool. 

The mapping technique déveloped in this thesis is based 
on the closed loop performance function, the poles of which 
are actually the roots of the system characteristic equation. 
The technique as developed is applicable to systems whose 
closed loop performance functions contain one and only one 
pair of complex conjugate poles. 

Now consider such a performance function consisting of 
one pair of complex poles and an assortment of real poles 
and zeros. The poles and zeros of the system can be plotted 
on a complex plane. If the locations of the real poles and 
zeros are fixed, a set of values for the transient response 
parameters defined in Chapter I can be determined for each 
possible position of the complex poles. If the values of 
each parameter are determined for a sufficient number of 


complex pole locations, the loci of constant values of each 


ange 





parameter can then be drawn on the complex plane. These loci 
are functions of the complex pole locations but are valid 
only for the fixed configuration of real poles and zeros. 

A simple example of the loci of constant values of a 
parameter exists in the case of a pure second order systen, 
that is, a system for which the closed loop performance 
function consists of ohne pair of complex poles and no real 
poles or zeros. In this case the overshoot is a function 
of damping ratio only and the lines of constant damping ratio 
are also lines of constant values of overshoot. These are 
radial lines emanating from the origin of the complex plane. 
This relationship is often displayed as a plot of overshoot 
versus damping ratio as is shown by Thaler (3). 

Clearly each configuration of real poles and zeros re- 
Sults in a different set of loci. For example, the addition 
of a real pole or zero to the second order case just con- 
Sidered results in lines of constant overshoot which are not 
straight lines. In fact, the exact location and shape of 
the lines of constant overshoot are determined by the actual 
locations of the real pole or zero. 

This leads to the need for a method of classifying 
systems according to the real pole and zero configuration. 
In order to facilitate the investigation of the sets of 
loci associated with different systems, the systems have 


been classified as follows: 
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a. Pure Second Order System - Closed loop 
performance function consists of two 
complex poles only. 


b. Second Order System With One Zero - 


Closed loop performance function con- 
Sists of one pair of complex poles and 
one real zero. 


c. Pure Third Order System - Closed loop 


performance function consists of one 
pair of complex poles and one real pole. 


These three classes of systems have been nondimen- 
Slonalized so that one set of loci of constant values of 
parameters is sufficient for all systems of that particular 
class. The details of the nondimensionalizing process and 
the complex plane plots of loci of constant parameter values 
are shown in Chapters III, IV, and V. Use of the families 
of loci determined in these chapters requires that the closed 
loop performance function of a system be known in factored 
form. 

Only the second quadrant of the complex plane is shown 


for two reasons. 
a. Only stable systems are considered so 


the poles are constrained to be in the 
second and third quadrants. 
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b. The complex poles are conjugates so 
the third quadrant loci are images of 
the second quadrant loci. 


In Chapter VI the technique is extended to third order 
systems with one or two zeros by means of an algebraic 
manipulation which permits evaluating the parameters of 
two or more systems of the classes listed above and com- 


bining the results. 


ile 





Chapter III 
THE PURE SECOND ORDER SYSTEM 


The second order system is the simplest type to which 
the previously outlined mapping technique is applicable. 
Its characteristics have been extensively studied and almost 
every servomechanism text contains a complete description of 
its transient response. The details of this transient re- 
Sponse will not be repeated here. 

In this chapter we are concerned with what may be called 
a pure second order system, that is, one whose closed loop 
performance function consists of only two complex conjugate 


poles. The performance function of such a system is 
2 2 
PR(s) = G48} = a 4h (1) 
(sta)“ +b 
The poles of this system are shown on the s-plane in 
Figure III-l. 
When R(t) is the unit step function the transient re- 
Sponse of the system becomes 


2.2 
C(t) =17+—18 tb 278 can (pt-o) (2) 


where O=arctan —— 


are 








FIGURE III-1] 


POLES OF THE PURE SECOND ORDER 
SYSTEM ON THE s-PLANE 


~18~- 





This is a very familiar expression, although it is usually 
written in polar coordinates, i.e., in terms of damping 
ratio and undamped natural frequency. 

It is convenient to define an artificial time unit 
as 


T = |jalt.. (3) 


The expression for the transient response in terms of T 


becomes 


e 
CUT) = 1 are on sin (SP - 9) (4) 


In terms of the artificial time, T, the entire 
transient response is determined if the ratio of a to b 


is known. This is equivalent to knowing the damping ratio. 


Time and Magnitude of Peak Overshoot 





The slope of the transient response is 


@ 2 2 

c(t) =S 42 e-8* sin pt (5) 
or 

e 2 2 

O(T) = at et sin on (6) 


From the nature of the transient response it can be 
seen that the peak overshoot occurs at the first time after 
time zero that the slope is zero. This requires that sin z 


be equal to zero, or 


=O 





(7) 


© is 
tt 
= 


(3) 


Therefore TT? es 


ol 


In real time 


C= (9) 


D 


1? 
i 
Oo} 


as is shown by Truxal (4). 

The loci of constant values of time of peak overshoot 
(T,,) are shown mapped onto the s-plane in Figure III-2. 
Although a scale is shown on this plot for convenience, the 
loci are independent of scale as they depend only on the 
ratio of a to b.. Hence any convenient scale may be sub- 
stituted for the scale shown. 

The magnitude of peak overshoot can be determined by 
substituting T, into the expression for C(T), Equation 4. 
The result of this substitution is 


- TB 
c(T,,) =lte (10) 


The magnitude of peak overshoot is the value of the ex- 
ponential term and once more dependent only on the ratio 


a 
b° ra 


M, = 100 e Pg (11) 


Loci of constant magnitude of peak overshoot are shown 


mapped onto the s-plane in Figure III-3. They display the 
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same information that is usually presented in the form of 
overshoot versus damping ratio plots, as has been 


mentioned. 


Rise Time 


Rise time has been defined as the first time after zero 
time when the value of the transient response is unity. In- 


spection of Equation 4 shows that this occurs when 


sin (21 ~ 90) =0 (12) 
or 
Er = c= 6 
br - 9 = (13) 
=~ % . 2 _b 
Therefore T, = 59 =p arctan — (14) 


Loci of constant values of artificial rise time (T.) 


are shown mapped onto the s-plane in Figure III-4. 


settling Time 


Settling time has been defined as the time when the 
height of the transient envelope becomes two percent of 


the steady state output value, This occurs when 


a2 
“2 rie = JOD (i'5)) 
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Therefore 


T= sim 902 + = in (25 + 1) (16 ) 
Si: 2 be 


Loci of constant values of artificial settling time 


=.) are shown mapped onto the s-plane in Figure III-5. 


Use of the Curves 


To use the curves shown in Figures III-2 through III-5, 
it is necessary to determine a and b. Effectively this 
means that the roots of the system characteristic equation 
must be known. The method used in obtaining these roots is 
immaterial. 

Once the roots are known some convenient arbitrary 
scale can be chosen for the loci plots. Since these are 
complex plane plots the real and imaginary scales must be 
identical. The point -a + jb can then be located on each 


loci plot and the values of M., T_, T 


D p ne and Ts read off 


directly. 

The three artificial time parameters can then be con- 
verted to real times by dividing each by the value of a. 
The transient response may then be rapidly sketched by 
marking off known points, remembering that: 


a. The slope is zero at t = 0 


and t = to: 
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b. The transient oscillating frequency 


is = cycles per second. 


ec. The sine function is zero at t = th. 
and increasing with time. 


Figure III-6 is an example of the procedure for the 
performance function shown. An exact curve is included 


for comparison. 


Conclusion 


It has been shown that the loci of constant values of 


he, T Ts» and My can be mapped directly onto the s-plane 


»? r? 
for the pure second order system. In addition, these four 
parameters and the transient oscillating frequency present 
enough information about the transient response of such a 


system to permit a reasonably accurate sketch to be made. 
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Chapter IV 
THE SECOND ORDER SYSTEM WITH ONE ZERO 


The addition of a real zero to a pure second order system 
results in a different type of second order system whose 


closed loop performance function is 


PF (s ) = aé + b© . oo (1) 


d (sta)~ +b 


The singularities of this function are shown On the s-plane 
in Figure IV-1. The constant which is usually interpreted 
as system gain reflects the assumption of unity feedback 
characteristics. 

As pointed out by Truxal (4), the Laplace transform of 
the response of such a system to a unit step input may be 


written as 


bf) =a a= + bo i emis (2) 
se oe a Ne A eee 
e (sta)" +b q (sta)~ +b 


This response is the sum of the unit step response of a pure 
second order system plus a portion of the impulse response 
of the same system. The addition of the impulse response 
makes the system relatively more oscillatory than would be 


the case without the zero. The magnitude of this effect 
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FIGURE IV-1 


POLES AND ZERO ON THE s-PLANE 


Second Order System with One Zero 
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depends on the position of the zero relative to the 
position of the complex poles. 

It is also useful to regard the addition of the zero 
as a consequence of the use of proportional plus derivative 
control in the forward path of a conventional unity feed- 
back servomechanism. Another equally plausible view is to 
consider this case as an approximation of a higher order 
system, that is, one in which the singularities tnat are 
far out in the left half-plane are neglected. 

The expression for the transient response of a system 


whose performance function is given by Equation 1 is 


c(t) = + lente” \ (a-a)? + © sin(bt +@-0) (3) 
b 


b 
where @ = arctan Gzg? and 9 = arctan — 


The slope of this transient response is 


C(t) = ar tee \(a-a)® to© ee @ sin(pt +B) (4) 


We may now express a and b in terms of d by defining 


m, = 7a, (5) 
5c 2. (6) 


Figure IV-2 illustrates these quantities. 


Este 





FIGURE IV-2 
ILLUSTRATION OF THE QUANTITIES m5 AND m,, 


Second Order System with One Zero 
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Substituting these relationships into the expressions 


for transient response and slope of transient response, 


\ 2 2 2 2 
Ps m,~ +m, (1-m, ) + Mm, <<a 
2 


we have 








C(t) = sin (m,dt+f-9) (7) 
* a m, © + m, 2 a -m, dt 
Cie) = it, \(a-m,) +m,” e sin (m,dt +P) (8 ) 
m 
where 9 = arctan T-m, 
m 
a 2 
9 = arctan —= 


We now define an artificial time unit for this type of 


system as 


T= |alt 


In terms of the artificial time, T, the expressions for 


transient response and slope of transient response become 


2 2 2 2 
\m +m 1-m,)“ +m -m.,T 
c(T) =1 + 1 2 \\ = eg 1 Sin (m,T+2-9) (9) 
2 


2 2 
° m,~ +m,~) eae me T | 
C(T) = a (1-m, )* + m,” eit sin(m,7+p) (10) 


where 9 and @ are defined as before. 
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It is important that the reader visualize what has been 
done. The simplest view is that the transient response is 
expressed as a function of an arbitrary time, T, and that 
this response depends only on the position of the complex 
poles relative to the zero. One merely has to adopt the 
magnitude of d as the unit of measure and determine Mm. and 


m. by measuring a and b on the s-plane using this new unit. 


2 
In terms of T the transient response is then independent of 
the actual numbers associated with a, b and d. All that 


matters is the magnitudes of m, and m An example will 


a. 
illustrate the point. 


Consider the two closed loop performance functions 


pR(s), = Ole ft) (11) 


(s + .6)° + (.8) 


_  100(8'%+'100) 
PF (s ),, = a (12) 


(s + 60) 


In both cases m, = 0.6 and m, = 0.8, and in terms of T the 


responses are identical. 


6 


C(T) = 1 +12.118 eo" sin(.87 + 63.4° - 126.85°) (13) 


This response is shown in Figure IV-3 with three scales 

on the abcissa. The first scale is in units of T. The second 
is in real time and is applicable to PF(s),; where T = lt, 

or t,=T. The third scale is also in real time and 18 appli- 
cable to PF(s),, where T = 100t or t,= .O1T. 
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An alternate point of view is that the general per- 
formance function for this class of systems has been 


transformed to a new complex plane by the transformation 


Z2=m, + jm, = Tay (14) 
The purpose of the transformation 1s to place the zero 
at minus one on the z-plane. All possible second order 
systems with one zero can be transformed in this fashion. 
If the variable z is treated as a Laplace variable, use 
of the inversion integral yields an expression for transient 
response, not in real time, but in an artificial time, 

T = dt. The result of this procedure is identical with the 
result obtained by manipulation of the expression for C(t). 
As has been indicated, the transient response as a 
function of T is dependent only on the values of Mm and Mm,» 


Therefore loci of constant values of magnitude of peak 


overshoot and the time parameters T_, T 


p oe and Ts can be 


drawn on the z-plane. 

For convenience the z-plane may be regarded as the 
S-plane rescaled so that the zero is placed at minus one. 
This places the complex conjugate poles at -m, + jm,. 
Application of the Laplace inversion integral yields an 
expression in T which can then be easily converted to real 


time, 





Magnitude and Time of Peak Overshoot 
The peak overshoot occurs at the first time after zero 
time that the slope of the transient response is zero. 


This requires that 


sin (mT + @) =0 (15) 
or 
m 
T= = = gy - arctan y=) (16) 


Loci of constant values of artificial time of peak over- 
shoot (T,,) are shown on the z-plane in Figures IV-4 and 
iy -5. 

To find the value of the systemoutput at the time of 
peak overshoot we substitute the expression for T. into 


p 
the expression for transient response, Equation 9. 


= gone 
m, + m,* (1-m, )* + m,” . n> B) 


i sin(r-9) (17) 


c(t.) = 1] + 


This reduces to 


mm 
2 2 2 2 - —(r-9) 
m,~ +m 1-m,)“ +m m 
Vm? + 05° \o-m)? +m, So °” iho! 
2] 


=] + 
c(T,,) 
Mo 
Since, as previously noted, tan 9 = =a,” we also know that 
=o 
sin 9 = ——\—___ (19) 
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This can then be used to reduce Equation (18) to 
m 


= * a(t) 


¢ (T.,) = 1 + \\(a-m,)? +m,” e (20) 


The exponential term of this expression is the amount of 
peak overshoot. Magnitude of peak overshoot expressed as 


a percentage of steady state output then becomes 


- lig) 
M. = 100 (ian: )° me ne % (21) 
p : mi me Of 


Loci of constant values of peak overshoot (M,) are 
shown on the z-plane in Figures IV-6 and IV-7. 
It is of interest to note that Equation 21 contains 
the expression for pork overshoot of a pure second order 
1 a 


System. The ratio nm. is identical to the ratio - for any 


given values of a and b, so 


e7ma/bd x pO Be (22) 


It appears that the peak overshoot of a second order system 
with a zero is the peak overshoot associated with a pure 
second order system with the same complex roots, modified 
by a factor which depends on the location of the zero re- 


lative to these roots. 
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Rise Time 


The system output is equal to its steady state output 
for the first time when the sine function in Equation 9 
1s equal to zero. The argument of the sine is always 
negative at zero time and increases with time. Therefore 
the first time after zero when the sine function is zero 


is when the argument is zero. That is 


mT +G-902=0 (23) 
or 
_ 6 40 
T., = — (24 ) 


Loci of constant values of artificial rise time 6) are 


shown plotted on the z-plane in Figures IV-8 and IV-9. 


Settling Time 
settling time is reached when the coefficient of the 
sine term in Equation 9 decays to .02 or two per cent of 


the value of the steady state output. This requires that 


\ m, 2 + my° \(a-m, )? + m,° -m,T 
1 = ,02 (25) 
2 
A simpler expression for Ts results if it is noted that 


C(O) = 0, as demanded by the initial conditions. In this 


case Equation 9 becomes 
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m,“ +m, \(a-m, + M, 


ee = 
a 


sin (g-9) (26 ) 


Now Equations 25 and 26 may be combined to yield 


~m,T 


- sin(@=0) = .02 (27) 


From this the expression for artificial settling time is 


a ae a in | -.02 sin(g-0) | (28 ) 


Loci of constant values of artificial settling time, 
(T.) are shown on the z-plane in Figures IV-10 and IV-11. 
The loci are omitted in the zone where peak overshoot is 
less than two per cent because settling time as defined 


here is meaningless in that area. 


Transient R@sponse Parameters 
Along the Z-Plane Axes 
The values of the transient response parameters as 


either m, or m approach zero are of some interest as 


2 
limiting cases. It is useful to visualize the angles © 


and @ because they appear so often in the expressions for 
the parameters. They are shown in Figure IV-1le. 


First, consider the parameters as m, approaches zero, 


1 


that 1s,along the m, axis or, alternately, on the s-plane 


2 
along the jw axis. In the limit we have 


Sng - 





geese gers yee gegen event estore eager eee Eee Eee Espen ee SP TeTE 
pects etp fete geriseoMecy sot eee ee Tg HUW UHeted ieee ota etsed 2 eesed epee eee eee see deserted tet 
EE ee 
detec te gesetceetctesct deters Gee sage reece gE AE Gece SEE EERE SESE EE EES SE EEE EE EE EEE 
Esse esepepesgemoeece ge teen trie tree pected sea a] aeesaeeaesebaegeaaeic LE VEU an eee iss acy 
ee eee ee eS eer SSeeP TAT TERR ect ee feeeTt LzdEs TEER CeCe PPE Eee eC 


FSS ES 


rial 


faletetet sol oalcfoloeleilalaetateleate lo foal iain (oobi eta ane mel ner ett aio eed | eters) 
BE Oe 





sad eevee espe fee ve Engen cE Leet gece eee gigi eae eevee ded] eee eH eedTSC eee ee eee TER 
EC EEE ER EEE ee ECE EEE CECE EEE EERE EERE EEE CORE HSE ee so 
4 oo 8 BUSES BESET FECES OSS Ee een anil een Pt te Ht CEE EEE EEE EH EEE EEE EET Hoeft tT 
ESSE copes Eeesgeceet ete CsCeCSSoesTectacefectesttf eco eee CO Po So pope eH iia ma tag iiiaiacbiis tare aig en 





re 
TT eee Fe 


4 Bet Bt tH 
alae a[eloala) alae eel slelaia)alalaiaeleiaialaleaea es ea 


pe eee Sass See Ee eee err ee ere ee ee eee eae ee er eee ee eee eee 
Pao re eeeeerrter eee eee ree Eee eee eee eee reece eee ee ee ce eee 
Pee ec eee er reece eee 
ESE ECCS oo eee eee ECE SCE EEE EE EEE oor 
Fo eet eR ee tr tr RHEE ESR EERE SES EH eer et eb er ee et ee 
PA pete CI CECE EE BE ECE 
| Het + trt4t fetes feayea npn eee) | oS ese erse lel olela) ht fags ce] eo] ele la | es) aloes) alalmpala) | apapa fe) el aay emia] sympa pm) ela pay lel apenas |e lala 4 fae fe 2H fed ee 
a 
ry 








FRET 
oe te bee ay fe 4 = tf 4 Ee 
—— rs r Lt 


on 
P| 
Cy 
a 





H LT) | 
PEPE epee CeCe ere eer eee ere ee eee eee eee eae eee eee 
cote Pet peeemeetreeet tttere FEE EERE EEE EEE EEE EEE ECE EEE COE EEE ECE IEEE ea ee 
> HS HCE Peer eee eee eee Cee rate i ce eee ere eee eee eee eee pee ee eee 
eeARepetedeeet sree Eee agaeeaeeggeee ARSE tz aaa dT Penge USP Go eve c Ee EeeaeEEd Te E 
ro tH t+ een CH ia ae +t th = pte a 


Pa So Cee SEE EEC ee 
oe BERT ER HERE EEE EESEe CES Eee Eee eee eet 


a Pt eh hs pg tp bet tty fd tT ed tn Ree os 
Sd eee Se Hicks Lace Rese ee gece Se a aes 
Sree an nu neem Trea ae ae ee REE ea Cate ae ea eettn eet en ee an oe gee ee eve cet ene "HE 

Pe po CO CORSE eee eePPEREEEr eee EEE Ee EEE EEE PPR EEE ee Ege eee 





Soeeeeeerbererirrr eee tee eet HOC QOS Sa250 0088 ee Eee bape) | | tS on 
Ganra js) 


rer yy yy + TT a a Ley 
BREE ete ee EA a EG EE SRS 


eee ete TT) (enpeslesh =.) al@reie pel p ala) 





PEE EEE EEE EEE SHEE EET Spite eet tir ttt ee Dy a 
LS fe ee Bsa Lake semeapeterete eect LUE accede dee eed LS 
eonea|cetegisceducredarqargatdun {iran ianjatiuedangeeiusais iuequaduequats cdeedeaqeeinse sails @isupen ial tsueatecued 
M SEES EEE EEE SS EEE EEE EEE EEE EEE EEE ET 
SESE ea ARE ERASER ESHER 
BE Eo SR EEE ERR SEE EEE PEE EEE COC 
SSIS SSS eee ee ater ce eect sen ed 
SSE eee eee ae au eee at ae ee eee 
BES E EEE EERE EERE EE EERE EERE CEE EE ECR EEE EE EERE EERE EEE ECE EEC EEC CCEE COE COS EEC EEE Ee 
Ea ar eee ete ere Me or Re er tte tt EE ae Meg Se 
EEE Eee CCN 
ESSSeESrSeeeeeE red eet fost cuE acta bepeemeereny erg eter eter eer FECES 
PT Se oon. PTT Sa i 9p fo eee [Fi H 


beer eeg reer cerry Cee CHEE Hee EEE EEE ES SEE EE PTET SEEN EEE EA a 


it | J cy B 
AAA a EAN HH pt alee) 
HEE EEC ERE 
| 
BEE COCCEeCCRO eal Cee 





TT Hy +4 YI iS 
aleleletaleya A HH tt tre fy 
fae] alngete |S] [aya ya ole § Heeegee 2 





subsnee  ¢auceslususGuguseuecs cues Poe PEER EEE EEE So HH PEE ECA EHEC SECS 
poeee REE Re 


Sete pasa Seieiers PeiietHee PSEC 





i EitbenaststasceceenintRatt 
ESeeeeceeee at aed td EE ECE ESE PE EERE 2d SEES ESS Seog ee 


SS Po Snes pao aos PSS [eee aE persis 





ERE Ee ELEC CEES EEE 





BATFE EEC EEE EEE REECE EEE EEE EEE EEE EEE EEE EEE EEE EEE 
eset ce Seg ee ese secs seats cte ecclissi SSH Ee agit taedneGEEE 
essere ensue bebe epg ae gued shes cae sataed fan irag eal suis seaee ecenae setae ust Gesaegeatdectesd sates 
S| ease 05 eee ee feohe: [eam See Reece Mes ree! Soi mene eles 8 Peers) eh ep ape pel eye) spe age Spee e Spa ele fe eee ee) ee ala ee ee aS 2) aie eee ete else eels oa le Sasa ab lees la ale pepelalofapelet=(- tsar. ]_] 
eset s oes ens ee Senge creed ae tesge gece ain GE ged ae gus Ege ees EEE ERE reg eee EEE 
FES SSE SHES ene ceeg eee gene eee a eee ee EE eee eee eee Seer peeeeteeeeeccet 
cope eee eee bres SE rr cept eer ee ee eee ee eee eo 
ere eed ee Sse Eae seg Egre gsc SE Gaee Ege Ces EEE SEER EGE REESE 
cesses sees eh eden guage age snes dE EGE ReGen ESE ELS EEE EEE Ene Sn ee EERE EEC 
ot tH can HH eves Vrystele eee ae eis aa) eae |) nooo eee Stata Benes Bee) ett tah halite! ptt Rt iH 














HHH fe cee fm ca af ae ffs fe fens) faa 

pee Be reece eet pp rr pm fp fp 
a Tittde tite eee Seed F 

FEST ferrprmerbersttces tte reet errr rere Ee etree ree Eee eee eee reer ee Pe eee eee eee eee ee Pe oe 

auers eS 

Gneia 

Ho 


BE Ce Eee eee ee eee eee eee ee eee eee ee ee eee ee eee eee eee 
SP ee SEE asad ia gs bese eee Ret! 
4 bo ra te bt ef ated bt ttt Pt ro Pi b+ b+ Ly +44 ts | Peete pases oo Lee bt + | +++ +4 44 aaeen po HEH} Seana HH Pet Pit 5c] a Rot Ho H+ Lt oH bot ire po 
eee ee 





eee eee ge 

ese eee eee gee eee 

SAS aecsenere reece geese ge Heer eee ese neg eg Rec gered EE ESE Sete 

cskcefesess aiff este ale aed dt eestee cree geese HERES EGG Gee gee se eS Ea 

Poe See reece epee ere Eee eee BEERS eee 
EEE 








Se eee 
Sa er fesse eee emer reer eee eee eere ee ee reer eee eee eee eee ee Ee EEE Eee eee eee 
re BL oe sheers PTT | peter ep oe pope ttle) a) Greys D6 Fe fon on pf op end ed at Hamp eign op caf nf caf Mapa fi ou] Plsep ra uf es) RPT G mH Ff ON capa] Ong YN ag emf Oped CNN smo UNL] MN ma LY waa WN Ye] OY ef cof ae] Fy cs] ey sp |] Pep alo] | feed oy oo] ee Cea ete |p| 

RC eee 
ee ee ee eee ease 
7 SST ep Te EE ep 
EH feeeteenqeseereesd spate tedeseetiget see icd ee agen ste ad speHAeUeh ed EES PE Hts diets eer ada Feet ee ePaper 























SPSS TESA SU 2S SEA SSE Skee Snag SSR eae eed ode ST 
SEES SSS SESS SESE a eer eae gee SSE ep Se deep eee edge 
Et ag fesecdaeegeere eee ered eee ecLeet gE EEE EL cne ga Et geet eee EEE GEE EEE nee eee EE 
Seed eed ceege ae ageie LGR HN aga GEES Obeeapee opp eErepTEEESE ce SUH PE SEET SU Ege REE EERE EEG Ee REE ECE eee eg 
SEE eer 14 FeeeEeT EL Sate eepep eH EEE EEE EEE EERE ERTS 
een eget] TERE EeeTEE PES Ered CePEEEEEEeEEeezE EEE EET EE EET EEC EEE EEE eee 
Feeder ere] eed Pe eeeen ay ceed Cea uot eetereng Trees Freezer EEE EE Eee EeE EEE EE EEE Lg EEE PEE EEE ETE Eee Pere 


Si ea pe nO rR Se er ree 


| SPS SS Seppe Ua peed it EC ESE SE 








= i aeeeeec reece Cee eeeeececet Etec Eee eee tet eee eee eee ELEC eeee eeeENEeEEeeEeEEEReeeREeeeESeeS Ree Ree eg RE 
ede PTE tle ea aa ada veneered ede eg eg Ee SE EE 
aehES 


EEE PEE ES eee ee eee eee ee Seerieettcieastbateitegs 


fay imygea ei) H+ Bee ete td F#RSR Osee8 
ECE Bo EH HEHEHE Ete Eee PEE z : 


SHEERS EES 


CECE ECE CECE EE 





ae plane 
“Mit jm. jie 


nae Mo 


FIGURE IV-12 
ANGLES 9 AND @ ON THE z-PLANE 
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on) eat 
Ty, = mi (rf - arctan m, ) (29) 
2 

M, = 100 \2 +m, % (30) 
T= 2 arctan m, + 4/2 (31) 

Mo 

Ts is infinite. 

As m., approaches zero there are two distinctly dif- 

ferent possibilities. For ms less than one, in the limit 
as m, approaches zero, Ty? Tw and Ts become infinite and 


M, becomes zero. 


On the other hand if Mm, is greater than one, in the 


limit as m,, approaches zero, we have 


*p ~~ I-my m, (32) 
af, 1-m, 7) 
. 100 (1-m, ) e 0 (33) 
i 
ae = = (34) 
1g my 7 


T. is infinite. 


These expressions are useful in determining the in- 


tercepts of a particular locus, 
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Construction of the Loci 


Each of the families of loci were constructed by cal- 
culating and plotting the values of the parameters as 
functions of Mm, for a number of fixed values of m, - The 
result was a series of cross-sections of the surfaces de- 
fined by the parameters over the z-plane. The coordinates 
my and m, were then read from these cross-section curves 
for each locus desired, and these points were plotted on 
the z-plane. This method was required because no method 
could be found to determine the equations of the loci on 


Use of the Curves 


The curves are used to determine the transient response 
parameters in the same manner as outlined in Chapter IIT 
in connection with pure second order systems. Once the 
closed loop performance function is known in factored 
form, that is, the quantities a, b, and d are determined, 
m and m., mast be calculated. This is relatively simple 


since 


a 


The point -m, + jm, is then located on each family of 


EB. . oe and M. are read off directly. 


loci and values of Tp? yr? tgp 0 


ao Sie 





The three time parameters can then be converted to real 


time by dividing each by the value of d, i.e., 


ao 


oe 2 S 
t= t=s— ay 


a ~ 
D ray’ _ dy’ and t. = 


With these three times and the magnitude of peak over- 
snoot the transient response may be rapidly sketched if it 


is remembered that: 


a, Tne slope at t = ty is zero, 


b. The transient oscillating 
frequency is or cycles per 
second. 


c. The sinusoid is zero at 
t= ty and is increasing 
with time. 


ad. The slope at zero time is 
BETS 

as which can be easily 

shown by use of the initial 


value theoren. 





Observations 


The loci of constant values of peak overshoot become 


Straight lines as they approach the origin. If m, and m,, 


are much less than one, the expression for peak overshoot 


degenerates to m, 


m 
M =100e © & 
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As has been shown, this is identical to the expression for 
peak overshoot for pure second order systems. These loci 


are actually very nearly straight lines if m and m, are 


2 
less than about 0.6, indicating that a zero will not in- 
validate a pure second order approximation unless its 
magnitude is less than about twice the system undamped 
natural frequency. For very low values of both m, and m,, 
the loci of constant values of parameters have been omitted 
because the pure second order curves are very nearly exact 

in this region. 

A comparison of the loci of constant values of peak 
overshoot shown in Figure IV-7 with the same loci for pure 
second order systems as shown in Figure III-3 will illustrate 
the previous observation. The comparison will also amplify 
the phrase "relatively more oscillatory than would be the 
case without the zero". For fixed complex root locations 
the addition of a zero can cause a large increase in the 
magnitude of peak overshoot. In fact, peak overshoot of 
greater than one hundred percent is a distinct possibility 
for such systems. It should be noted that the addition of 
a zero does not effect the transient oscillating frequency. 

The initial slope of the transient response is not zero 
es was the case with the pure second order systems. The 


initial conditions are met, however, because there is a 


velocity discontinuity at zero time due to the combination 
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of the input function and the zero. AS was mentioned, 
the zero can arise from the addition of derivative 
control in tne forward path. Differentiation of the 
unit step through this channel yields a unit impulse, 
and it is this impulse that accounts for the initial 


Slope. 
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Chapter V 
THE PURE THIRD ORDER SYSTEM 


The addition of a real root to the pure second order 
system results in a new system which we define as a pure 
third order system. The closed loop performance function 
of the pure third order system then has one pair of com- 
plex conjugate poles, one real pole and no zeros. 

Much less is known about the pure third order system 
than is the case with the pure second order system, This 
is due largely to the introduction of a third variable in 
the equations due to the third root. The mathematical 
expression for the transient response of the pure third 
order system contains a decaying exponential term in ad- 
dition to the constant and damped sinusoidal terms of the 
second order system. The actual expression for the 
transient response is not difficult to obtain, nor to plot, 
but rather the actual plotting is laborious. The main 
difficulty in establishing standard sets of curves for the 
third order system, as is often done with the. second order 
system, is in portraying these curves as functions of 
three variables, 

This does not imply that nothing has been done in 


this area. To the contrary, much effort has gone into the 
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study of third order systems. Unfortunately, the results 
of this effort has often produced new ways of factoring 
the system characteristic equations without producing any 
simple method of evaluating the transient response. One 
such method was developed by Y. J. Liu (5), and is based 
on nondimensionalizing the characteristic equation. By 
using this method, stability and the roots of the char- 
acteristic equation can be determined from a set of non- 
dimensionalized charts. 

Another method of importance in any discussion of 
third order systems is Mitrovic's method (6). In this 
method, by use of a selected transformation, it is 
possible to compute a set of curves which inherently con- 
tain the solution to all third order equations and can be 
used to determine stability. These curves have been com- 
puted and are known as third order charts, They may be 
used directly in compensation by cancellation provided 
suitable root locations are known. However, this still 
does not solve the problem of where to locate the roots 


and it does not yield transient response directly. 


Transient Response 


A plot of the roots of the pure third order system is 
shown in Figure V-1. The closed loop: performance function 


for the unity feedback third order system may be written as 
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FIGURE V-1 
POLES OF THE PURE THIRD 
ORDER SYSTEM ON THE s~-PLANE 
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2 2 
Eis) & a +b oe 
PF(s) = aol ste ) (statj sta-j (1a) 


This may be simplified as follows 


= a“ + b“)c 
men (ste) |(sta) Sel om) 


If a unit step input is applied to a system having 
the performance function shown in Equation 1, the complex 


frequency domain output becomes 
2 2 
ot = pl tt Pe (2) 
s(sic) [ (sta) + bd°| 


The inverse of this transform gives the time domain out- 
put of the system, which is 


c(t) =1 - AP, ost 


a-c)" +b 


ae 2. s 
4 Set be 68 san (pt-0-Z) (3) 


b |(c-a)? + b° 


where 


b 
9 = arctan > 


2 b 
§ = arctan > 


Now if we redefine the root locations in terms of the 


location of the real root as shown in Figure V-2, the real 


aoRe 
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FIGURE V-e2 
ILLUSTRATION OF QUANTITIES ny AND no 
PURE THIRD ORDER SYSTEM 
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and imaginary parts of the complex roots become 


a= n, |¢| 
b = n,|c| 


The quantities n, and n, are merely rescaled values of the 


1 2 
réal and imaginary parts of the complex roots. 


Making these substitutions for a and b in Equation 3, 


the expression for transient response becomes 


e 2 


n + Ne -ct 


fan <<) - — 


(1-n, ) + a, 


n tn -n,ct 
v A Z e t sin(n,ct-9-9) (4) 


\ 2 2 


9 = arctan np/-n, 


where 


Me 
@ = ama Ta 
This substitution has partially nondimensionalized 
the transient response of systems with the same relative 
root configuration, The transient response can now be 
completely nondimensionalized by making the additional sub- 
stitution T = |c|t. This is equivalent to defining a new 
artificial time for this class of systems as was done in 


Chapters III and IV. The transient response now becomes 


AG6< 





+n 2) 


2 
(ny 2 


-T 
(1-n, ) +n 7 5) 


c(T) =1- 
2 


1 2 1 sin (n,T-9-2) 


2 2 
nN, | a-n,) +, 


It is of interest to look at a different derivation 
of Equation 5. As a first step let us write the closed 
loop performance function in terms of the root locations 


shown in Figure V-2. 


2 aes 
PF(s) = ea (6a) 


ea Vealctt tne! 
(sic) [(sin,c)* + (n,c)* | 


Dividing numerator and denominator by o? 


gives 
rn e +n . 


(= +1) | +n,)* +n," | 


Now if we make the linear transformation 


ef ese 
+ jn, = 


\c| 


2 = ny 


the performance function becomes 


2 2 
no +n 


| 2 
PF (z) (ea) [(ein,)? +n," (7) 


This has effectively transformed the system onto @ 


new complex plane in which the real root is located at 
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minus one and the complex roots are at -n, + jn,- Now 
upon application of a unit step the complex frequency 


domain output becomes 


C0) ee (8) 
z (zt) | (24m, ) + | 


The inverse transform of this gives the transient response 


as 
2 2 
n.~ +n 
c(t) =1 - the 
1-n +n 
. (9) 
2 2 
n +n -n,T 
} ee eee 1 sin (n,T-9-9) 


2 2 
Ne \(2-n, ) + ny 


In this case T is the time associated with the inverse 
from the z-plane. Note that Equation 9 is identical to 


Equation 5. 


Time and Magnitude of Peak Overshoot 


As in the previous cases, the slope of the transient 


response is zero at the time of peak overshoot. That is 


c- C(t.) = 0 (10) 


The Laplace transform of this for zero initial conditions 


is 
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> ES o(t)| = S C(s) (11) 


Multiplying both sides of Equation 2 by s gives 


S C(s) eg ha (12) 


(ste ) [ (s+a) + b 


Combining Equations 11 and 12, and obtaining the inverse 


transform gives the slope of the transient response as 


2 2 -ct 
d = (a tebmje p 
ae CMe) fa tee 


(c-a)~ +b 








2 2 ~at 
+ (24D) _e P sin(vt,-4) = 0 (13) 


b \(c-a)? ne 


where @ is defined as before @® = arctan a 


Equation 13 can be reduced to 


-ct Vea )@ane -at 
c (ac+ b@)lb e€ ee (c-a)=+b° e Pgin (>% = g) | = 4 (14) 


b | (c-a)® + Det 


For practical systems the left side of Equation 14 can 
equal zero only when the term in brackets equals zero. 
Therefore, the condition to be satisfied at time of peak 


overshoot is 
eu, -at 
be (th)-an-e 2 sin (bt, - 6) = 0 (15) 


-63- 





Equation 15 can now be nondimensionalized by making the 


same substitutions as used in the transient response, i.e., 


a=n,|c| 
b = n,| ¢ | 


< 


le|t 


When these substitutions are made in Equation 15 the con- 


dition to be satisfied at time of peak overshoot becomes 


——s pia = sin(n,T - @) (16) 


| (2-n, )* +n, 


"2 
where @ = arctan I-n, 





From Equation 5 we see that the magnitude of peak 


overshoot (expressed as a percentage) can be written as 


2 


2 
n +n -n.T 
M, = 100 1 € e iP sin(n,T ,-0-2) 
2 2 
n, | (1-n,) +n, 
2 2 r (i) 
ny + Ns 2 P 7 
oO 
wa “ny an no" 


For any combination of n, and n,; Equation 16 can be 
solved by trial and error techniques to determine the 


values of T for the first, second and each succeeding maxi- 


Wale 


mum of the transient response. With a small amount of 
practice these values of T can be determined with one to 
three iterations per value. The values of Tp determined 
from Equation 16 can then be inserted in Equation 17 to 
determine which maximum gives the peak overshoot, and 
the amount of that overshoot. From the foregoing we see 
that the calculations for time and magnitude of peak over- 
shoot are inseparable. Each must be known before the 
other can be determined correctly. This is not a severe 
handicap, however, as the peak overshoot occurs on the 
first maximum for practically all combinations of n. and 


1 
n. which include values of n. less than two. 


2 2 
Loci of constant values of magnitude of peak overshoot 
(M,) are shown on the z-plane in Figure V-3. Another set 
of loci of constant values of peak overshoot with a com- 
pressed ordinate scale to permit extension of tne plot to 


include larger values of n, is shown in Figure V-4. Although 


2 
Figure V-4 does not represent a complex plane due to the 
adifference in scales used on the axes, it may be considered 
as one for the purpose of reading values of peak overshoot 
associated with values of ny and Np: 
Loci of constant values of artificial time of peak 
overshoot (T) are shown on the z-plane in Figure V-5. 
2 
increases the peak overshoot occurs on the second or sub- 


An interesting point in connection with Th is that as n 


sequent maximum of the transient response curve, and the 
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time of peak overshoot is no longer unique. The zone in 
which this situation exists is shown in Figure V-6 which 


is an extension of Figure V-5. 


Rise Time 


Rise time is the first time at which the output of a 
system equals its steady state value. This is equivalent 
to saying that rise time is the first time for which the 
overshoot is zero. From Equation 5 we see that the con- 


dition for zero overshoot is 


¥ + ny” _T n,* os no -nyT 
eT = +S e + sin(n,T-0-f) (18) 


fom +n. \(.- 2 2 


This can be rearranged as follows 


n 


sin(n,T-0-f) =n (19) 





Obviously, this too is a transcendental equation, but 
just as in the case of time of peak overshoot, it can be 
solved by trial and error methods for any selected com- 
bination of ns and No > 

Loci of constant values of artificial rise time (T..) 


are shown on the z-plane in Figure V-7. 
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settling Time 


As normally defined, settling time is the last time at 
which the output error exceeds some fixed percentage of 
the steady state output. This error is often specified as 
two percent. Using this definition, settling time is ex- 
tremely difficult to determine for a third order system as 


the condition that must exist at T, is 


2 2 2 2 
n + mn -n.T +n 
ae ee e : sin(n,T-9-9) ~ —+ £5 et a + 02 
ny (1-n, )* a — (1-n,)" +n, 


(20) 


To determine settling time exactly would require several 
trial and error solutions of Equation 20 or a plot of 
Equation 5 for each combination of n, and ni, The labor 
involved in either method is prohibitive. 

The settling time for the third order system as defined 
nere is the time at which the sum of (1 - exponential term) 
plus the height of the envelope of the sinusoidal term is 
equal to two percent of the steady state output. The 
equation to be satisfied for this condition is 


2 2 ) 2 2 
-n.T n +N, -T. 


n +n 
—————- "a = .02 (21) 
l-n +n 
n, \@-n,) + Ny 1 2 


Loci of constant values of artificial settling time 
(T.) are shown on the z-plane in Figure V-8. Settling time 
| = 
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for systems with less than two percent overshoot is not 
shown because settling time as defined here is meaningless 


for such systems. 


Construction of the Loci 


In constructing the loci of constant parameter values 
it was necessary to calculate the values of each parameter 
for a large number of combinations of n, and n,- This was 
done by calculating the parameter values as a function of 
no for a number of fixed values of nj: All of the calcu- 
lations were performed by trial and error as previously 
indicated. 

The values of the parameters were then plotted as a 
function of n, or n, with the other variable n, or n, as 
a parameter of the curves. This gave a series of cross- 
sections of the surfaces defined by the transient response 
parameters over the z-plane. The coordinates ny and nN, 
were then read from these cross-section curves for each 
locus desired. These points were plotted on the z-plane 


and the constant value loci drawn. 


Observations 


The loci of constant values of peak overshoot become 
straight lines as they approach the origin of the complex 


plane. This is consistent with the usual second order 
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approximation for the case in which the value of the real 
pole is several times the value of the undamped natural 
frequency of the system. In Figure V-9 the validity of 
this approximation is demonstrated by a comparison of the 
peak overshoot of the pure second order system with that 
of the pure third order system for small values of n, and 
n° 

Another area of interest on the plot of overshoot loci 
is the area of zero overshoot. The zero overshoot locus 
represents the boundary where the magnitude of peak over- 
shoot becomes too small to calculate. The boundary as 
shown represents a value of the order of e719. For values 
of n, greater than about 0.94 the peak overshoot is too 
small to calculate, regardless of the value of no° In this 
case the response is essentially that of a first order 
system with a decaying ripple superimposed on the rising 
exponential curve. 

As noted previously, there exists a region in which 
the time of peak overshoot is not necessarily unique. 
This area is shown in Figure V-6. The zone is bounded by 
a line through the values of ny and n,, for which the 
height of the first and second maximums of the transient 
response curve are identical. Obviously this includes 
the entire imaginary axis. 

For small values of No the time of peak overshoot is 


relatively independent of the value of ny: In fact, the 
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only area where time of peak overshoot becomes sensitive 


to the value of n, is where both n, and n, are large. 


1 
For small values of n 


2 


2 the rise time is also relatively 
independent of the value of n> although it is more dependent 
om this value than the time of peak overshoot is. For large 
values of both ny and no the rise time becomes more depend- 
ent on the value of n, than on the value of n,- 


For small values of n, the settling time is essentially 


iu 
independent of the value of No» and settling time is never 
as dependent on the value of n, as it is on the value of ny ° 

To avoid heavy concentrations of lines the loci are not 
plotted near the origin of the complex plane. This does 
not detract from their usefulness as the second order 


approximation is valid in this area. 


Use of the Curves 


The curves for the pure third order system are used in 
exactly the same manner as the curves of Chapters III and 
IV. Since their use has already been explained there, 
the explanation will not be repeated here. Instead an 
actual example is presented. 

Assume that a system has the following closed loop 


performance function 


2500 
PF (8) = (SaEy(EeSeIE) HOTTY 
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Dividing the real and imaginary parts of the complex poles 
by the value of the real pole gives: 


= 


=] 


UY Uiln 


= 8 


S 
Mm 


Entering Figures V-3, V-5, V-7, and V-8 with these values 


of ny and n, we can read the following parameter values: 
M = 13% 
vl = Syl 


B.C 


r 
tt 


T. = 10.2 


Dividing the three artificial times by the value of the 


real pole gives the real times: 


_ 521 7 
t van 5 . 1,02 8EéC,. 
S40 
Ben = = ,76 sec. 
_16.2 << 


From the performance function we see that the transient 
oscillating frequency is four radians per second, or the 
period is = seconds, The three values of time, the period 
of oscillations, magnitude of peak overshoot and the fact 
that the slope is zero at zero time and at time of peak 


overshoot constitute sufficient information from which to 
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sketch tne transient response. Tne transient response is 
sketcned in Figure V-10. 

Now if we apply a unit step input to the system and 
obtain tne transient response by performing the inverse 


transformation, we get 


-5¢ 


C(t) = 1 - Be" + 1,11607@*sin(229¢ - 169.7) 


This output has been calculated and is plotted as a function 
of time in Figure V-1l. From an expanded scale plot of this 
output, transient response parameters have been determined. 


These are compared with the predicted values below. 


Predicted Value from 
Parameter Value Transient Response 
ty 1.02 sec. 1.025 sec. 
te. .76 sec, .76 sec. 
t. 2.04 sec, 2,00 sec. 
M, 14 &% 12.7 &@ 
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Chapter VI 
EXTENSION OF THE TECHNIQUE 


The curves previously developed can be used to approxi- 
mate the transient response of third order systems with 
one or two zeros. This extension is made possible by a 
Simple algebraic manipulation whicn breaks a system down 
into two or three parallel systems. The transient response 
of each of the parallel systems can be approximated using 
the curves of Chapters III, IV and V. These approximate 
transient responses can then be combined to give an approxi- 
mation of the transient response of the original system. 

The scheme can be explained best by an illustration. 


Consider the closed loop performance function 


PF(s) ~ ola bh). —__test) (1) 


(st+c) | (sta)~ + b 


This can be rewritten as 


iii, S(t) Stet? (2) 


q y (s+c ) | (sta) +b | 


where f is defined by 


f= (d - c) (3) 
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We can then write the performance function as 


2 2 2 2 
PF : a+b f (a~ +b . 4 
eo) d (s+a)* + “2 "a (stc)[ (sta) +b-] 7 


Applying a unit step input and taking the inverse transform 


gives tne transient response as 


C p= second | f ae third | 


C(t) = & | order response * @ |order response 


To sketch this transient response the pure second order 
response and the pure third order response must both be 
sketched first by the methods previously described. Both 
must be sketched as functions of real time. At any instant 
the total response is 5 multiplied by the value of the pure 
second order response at that time plus = multiplied by the 
value of the pure third order response at the same time. 
The two responses can be combined most rapidly by graphical 


methods. 


As an example, consider the performance function 


- 26 (st2) 
~~ (st) [(s+2)° + 3°| = 


In this case f = (d -c) = (2 - 4) = -2. Therefore 


a - 2 52 6 
PF(s) = 3 eae oe 2 aera (6) 


c(t) = 2 | Pure second | _, | Pure third 
order response order response 
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For the pure second order response, a = 2 and b = 35. 


From the curves in Chapter III values of parameters are: 


Dake ool 
T = 1.47 
T. = 4,1 


M = 12.8% 
Dividing the artificial times by the value of a ylelds: 


t. = 1.05 sec. 
p 5 


t. = .735 sec. 


Ce =m 2.05 sec. 


For the third order response: 


RT = 5 a +2 
™ = = = 


From the curves in Chapter V values of parameters are: 


56 


J 
H 


T = 4,3 


tj 
Hi 


8.5 
M = 6.5% 
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Dividing these artificial times by the value of c yields: 


f Sieh 4 sec, 
p 


cA =z 1.075 sec. 


t. = 2.125 sec. 


Note that in both cases the transient oscillating frequency 
1s = cycles per second. 
The process of combining the two responses is accomplished 
graphically. This is illustrated in Figure VI-1 where the 
two transient response curves and their sum have been sketched. 
The same method can be used to approximate the response 
of third order systems with two zeros. This merely involves 
more labor and of course the accuracy declines with the ad- 
Gition of zeros. 
Consider the following closed loop performance function 


Ga De) (std, ) (std, ) 


ete) = TS, Tore) [lead + OT 7) 


Tnis can be rewritten as 


2 2 s+c+f, ) (s+c+f 
PF(s) = et?) Bike Aes (8 ) 
1~2 (s+c ) |(s+a) + b | 
where 
fs as qd, - Cc 
f., = d., - Cc 
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Multiplying the terms in the numerator yields 


2 
mark o(at + 0°), (stc)~ + aia +f, f, (9) 
Pe (ste ) | (s+a) + b | 


This can be broken into three terms as follows: 


ES) wale: bo S+c a ffi tfel® | ae? ay® 
: 2 dd 


jd 190 (s4a)©4p* 


rey Baye Tata Ba 
(10) 





é 


1-2 (stc)|(sta)~ +b 


is fits c(a* + b°) 
182" (are) [(eva)™ +] 


This can be sketched in the same manner as before except 
that the transient response curve will be the sum of 
three curves instead of two. 

No satisfactory method of approximating the transient 
response of more complicated systems has been found. The 


only alternatives seem to be: 


a. Approximate higher order systems 
with lower order systems by 
ignoring the poles and zeros that 
are far out in the left hand plane, 
or 

b. Expand the performance function by 
the method of partial fractions into 
fragments of manageable size, utilizing 
the techniques developed in this thesis 
whenever possible. 
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Chapter VII 
CONCLUSIONS AND RECOMMENDATIONS 


The servomechanisms investigated in this thesis have 
been divided into classes according to the number and type 
of poles and zeros of their closed loop performance 
functions. This was done to facilitate the investigation 
of the transient responses: of' the different classes of 
systems. Systems investigated in this thesis are those 
whose closed loop performaance functions contain the 


following pole-zero configurations: 


a. Two complex poles only. 
(Pure second order systems ) 


b. Two complex poles and one 
zero. (Second order system 
with one zero) 


c. Two complex poles and one 
real pole. (Pure third 
order system) 


It has been shown that for each of these classes of 
systems the transient response of systems with the same 
relative pole-zero configuration can be nondimensionalized. 
The transient response can be expressed as a function of 


an artificial time and two complex coordinates. For the 
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pure second order system the artificial time is the magni- 
tude of the real part of the complex poles multiplied by 
real time and the complex coordinates are the real and 
imaginary parts of the complex poles. For systems with 
three singularities the artificial time is the magnitude 
of the real pole or zero multiplied by real time and the 
complex coordinates are the real and imaginary parts of 
the complex poles divided by the magnitude of the real 
pole or zero. 

With the transient response expressed in this manner, 
magnitude and time parameters of the response can be 
expressed in terms of percentage of steady state output 
and artificial time respectively. Loci of constant values 
of four transient response parameters for each of the 
three classes of systems listed above have been mapped onto 
a complex plane. These parameters are magnitude and arti- 
ficial time of peak overshoot, artificial rise time and 
artificial settling time. For the pure second order system 
the loci are mapped directly onto the s-plane, while for 
the systems with three singularities the loci are mapped 
onto an auxiliary complex plane. 

The approximate value of these four parameters for 
practically any specific system of the three classes dis- 
cussed can be determined from these families of loci of 


constant parameter values. A rather accurate transient 
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response curve can then be sketched using these four para- 
meters, transient oscillating frequency, initial slope and 
Slope at time of peak overshoot. This has been demonstrated 
by example. 

Although the discussion in this thesis has been limited 
to systems which exhibit unity feedback characteristics, 
this 18 not a necessary restriction. The closed loop per- 
formance function of a non-unity feedback system can be 
written as a constant multiplied by the performance function 
of a unity feedback system. Regardless of the value of this 
constant, it will appear in the transient response as a 
multiplier of all terms. Therefore, the value of the con- 
stant will not effect the value of any of the transient 
response parameters determined in this tnesis. The transient 
response of such a system will merely have a steady state 
value different from unity. 

Although loci of only four transient response parameters 
have been mapped onto complex planes in this thesis, loci 
of other parameters could be mapped in a similar manner. 
Other such parameters which would appear to be useful in 


improving the resulting transient response sketches include: 


a. Delay Time (Time at which the output 
first equals one-half the steady 
state value). 





b. Time and Magnitude of First 
Undershoot. 


Inclusion of additional singularities in the closed 
loop performance funtion of a system presents additional 
problems. Although the technique developed in this thesis 
can be applied to such systems, one set of loci is not 
sufficient for all systems of a class. Consider a system 
with two complex and two real singularities. Application 
of the mapping technique would require fixing one of the 
real variables as a function of the other real variable 
and construction of a set of loci for this fixed relative 
configuration of the two variables. Additional families 
of loci would be required for each fixed configuration 
of these two real variables for which knowledge of the 
parameters was desired. 

One method of avoiding the additional mapping required 
for third order systems with one or two zeros has been 
presented. This involves separation of such systems into 
parallel components whose individual responses can be 
approximated using the parameter loci presented in this 
thesis. Once the component responses have been approxi- 
mated they may be combined to give a resultant approxi- 
mation of the system transient response. Using this 
scheme the accuracy of the resultant approximation is 


considerably less than the accuracy of the component 
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approximations. However, careful construction of the com- 
ponent approximations should permit the construction of a 
useful resultant approximation, at least for the third 
order system with one zero. 

The technique developed in this thesis is not restricted 
to the mapping of transient response parameters, Frequency 
response parameters could also be mapped onto a complex 
plané by a similar technique. 

The authors acknowledge that they have not exhausted 
the possibilities of the technique presented in this 
thesis. In fact, they believe that there are a vast number 
of possibilities yet to be investigated. However, due to 
the labor involved in the mapping technique, further investi- 
gations should be conducted with the aid of a digital com- 
puter. Specific areas recommended for investigation using 


a computer are: 


a. The three classes shown in this 
thesis to extend the loci plots 
to include larger values of the 
complex coordinates. 


b. The third order system with one 
zero to include construction of 
sets of parameter loci for each 
of a number of fixed relative 
configurations of the real pole 
and zero. 
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Construction of loci of ad- 
ditional parameters such as 
delay time and time and 
magnitude of the first under- 
Shoot for the systems listed 
in a and b. 


Construction of loci of 
constant values of frequency 
response parameters. 


The possibility of mapping 
statistical quantities such 
as integral error squared 
for various inputs. 
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